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Some cosmological consequences of first order quantum corrections to Maxwell electrodynamics 
are investigated in the context of a spatially flat homogeneous and isotropic universe driven by a 
magnetic field plus a cosmological term A. The introduction of these quantum corrections may 
provide a more realistic model for the universe evolution. For a vanishing A, we derive the general 



solution corresponding to the particular one recently found by Novello et al. |gr-qc/9806076|. We 

also find a general solution for the case when A is a non- vanishing constant. Both solutions describe a 

^\. . non-singular, bouncing universe that begins arbitrarilly large, contracts to a minimum non-zero size 

0\ ' CLmin and expands thereafter. However, we show that the first order correction to the electromagnetic 

^\ , Lagrangean density, in which the analysis is based, fails to describe the dynamics near amin, since, 

7—i • at this point, the magnetic fields grows beyond the maximum strength allowed by the approximation 

<-». ' used {B <^ 8.6 x 10~^ Tesla = 0.0086 Gauss). The time range where the first order approximation 

flj , can be used is explicitly evaluated. These problems my be circumvented through the use of higher 



[j, • order terms in the effective Lagrangean, as numerical calculations performed by Novello et al. |gr 

qc/980908Cj, for the vanishing A case, have indicated. They could also be evaded in some models 
I~^ based on oscillatory behaviour of the fundamental constants. A third general solution corresponding 

to a constant magnetic field sustained by a time dependent A is derived. The temporal behaviour 
of A is univocally determined. This latter solution is capable of describing the whole cosmic history 
^ . and describes a universe that, although with vanishing curvature {K = 0), has a scale factor that 

(^ ' approaches zero asymptotically in the far past, reaches a maximum and then contracts back to 

OO . sm arbitrarilly small size. The cosmological term decays during the initial expansion phase and 

CO ' increases during the late contraction phase, so as to keep B constant throughout. An important 

Cn ' feature of this model is that it presents an inflationary dynamics except in a very short period of 

^^ , time near its point of maximum size. 
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I. INTRODUCTION 



C^ . In a recent paper, Novello and colaborators pi have analysed the cosmological consequences of quantum corrections 

L^' ■ to Maxwell electrodynamics. They have considered the quantum effects leading to the production of electron-positron 

pairs that were first derived by Heisenberg and Euler 0. The analysis was a semiclassical one (quantum field in 

^^ ' a classical general relativistic geometry), made in first order on the effective Lagrangean density (weak field limit) 

H ] and applied to a spatially flat Friedmann-Robertson- Walker (FRW) model. The more interesting result found in Rcf. 

■ - - ' was the removal of the primordial singularity due to the appearance of a negative pressure in the early stages 

of the universe. The analytical solution derived shows that the energy density associated with the electromagnetic 

field vanishes at the point where the scale factor reaches its minimum. The non-singular behaviour of the model is 

unaffected by the presence of other types of ultrarelativistic matter obeying the equation of state P(ur) = Piur)/'^- I^i 

a subsequent paper [pj, the analysis was extended beyond the first order approximation and the numerical solutions 

obtained show that the non-singular behaviour is preserved. 

The conclusion reached by the authors of references ||l|,p| was that the cosmological singularity of FRW models is 
a distinguished feature of classical electrodynamics. This problem is overcome when quantum corrections are taken 
into account, leading to a more realistic description of the universe. 

Lr the present paper we extend the analysis developed in [Q. We begin by showing that the analytical solution 
obtained in that paper is a particular one and write down the corresponding general solution. We then analyse how 
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this solution is modified by the presence of a non-zero cosmological term A in two different situations: we consider both 
the case of a constant A, which leads also to a non-singular universe, and the case of a time varying "cosmological 
constant" that supplies energy to a constant cosmological magnetic field. The equation that leads to a constant 
magnetic field as a possible solution to the field equations appeared already in Ref. j^] but was disregarded. The 
presence of a time varying A gives physical meaning to such a solution. The time dependence for A is univocally 
determined and represents a slightly modification of a form that has been suggested frequently in the literature 
(A = 37f^ — Co, Co = constant Q). The universe described by this latter solution begins arbitrarily small in the far 
past, expands towards a maximum size amax at the point where the vacuum energy reaches its minimum and then 
contracts towards an arbitrarily small size in the distant future. A remarkable feature of this model is that it has an 
inflationary dynamics, except during a short period of time near amax- This scenario obviously violate the standard 
evolution law for the magnetic field {B{t) ex a~^), which is commonly used to rescale to present time the several 
constraints on primordial magnetic fields presented in the literature |qM| . 

We also discuss the restrictions on the domain of validity of the solutions found, due to the use of the first order 
approximation for the effective electromagnetic Lagrangean. 

The paper is organized as follows: In Section |l| we set down the basic equations. In Section III we generalize the 
solution derived in lit] which assumes a vanishing cosmological term and a time-dependent magnetic field. . In Section 
IV we obtain a new solution that takes into account the presence of a constant A. A new solution involving the 
presence of a A(t) plus a constant magnetic field is derived is Section M. Section VI contains a summary of our results 



and suggestions for future work. 

In what follows Greek indices run from to 4 and Latin indices run from 1 to 3. As in pi], we use Heaviside-Lorentz 
electromagnetic units, but, unlike that reference, and unless otherwise stated, we make ?i = c = 1, . This is the 
scheme used in ^ which seems to be more appropriate to our needs. In this system, the magnetic field S = J? ]|^ is 
measured in Tesla |7| and the fine structure constant is given by jQ 

e2 1 

(1) 



47r 137.036 ' 
so that the magnitude of the electron charge is 



e = V4^ = 0.30282 . (2) 

(See the Appendix A of Ref. ]|j and the Appendix on Units and Dimensions of Rcf. Jgj for further details.) 

II. FUNDAMENTAL EQUATIONS 

In the system of units used in this paper, the Lagrangean density for free fields in classical Maxwell electrodynamics 
is written as 

^(MAXWELL) = ^-tF'^'^F^i, = — T-F , (3) 

where F^^ is the electromagnetic field strength tensor. Canonical energy-momentum tensor is then given by 

rpiMAXWELL) _ p ™ + ip„ (A\ 

-^ fill ^ /ja^ "^ zL '^'^ ' ^ ' 

The quantum corrections to classical electrodynamics that should be considered in order to describe the process of 
creation of electron-positron pairs by the electromagnetic field were evaluated by Heisenberg and Euler ]||,g] . A clear 
derivation of the effective Lagrangean that describes these quantum effects can be found in ]jiq] , where one should be 
aware that the Gaussian system of electromagnetic units is used. The first order calculation yields for the effective 
Lagrangean density p],p^ 



c = -\F + ^, 



}.F^ + —{F*f 
4 16^ ' 



(5) 



where 

F* = F*^Ff"' , (6) 

F*j^ is the dual of F^^, 



/^ 



" -■ 1.4x lO^GcV"" 



45 mi 



(7) 



and rUe = 5.110 x 10"'* GeV is the electron mass. [If one restore all units a factor h^ /c^ appears in (|^).] 

The domain of validy of the above expression for C is that of small wave numbers and small electromagnetic field 
streuffth, i. e., MM 



k <^me , 



(8) 



E <C Ecr = — 



2 ^2 
""e 



e ^/A^^a ' 



(9) 



B<^Bc 



e 



V4 7r( 



8.6 X lO"'^ Tesla = 0.0086 Gauss 



(10) 



[Restoring all units, a factor of c/Ti appears multipling the right hand side of (ph and a factor of c? /Ti multiplies the 
right hand side of equations (||) and (1^).] 

The corresponding modified energy-momentum tensor becomes [Oj 



dC 
T - -4-—F °'F 



( dC 



'"' dF ^ """ ' \dF'' 



F* -C 



9t^u 



We will aply the above equations for a homogeneous and isotropic universe with line element 



ds^ ^dt^ -a\t) 



1- Kr^ 



r^d0^ +r^sm'^i 



(11) 



(12) 



where K = 0,±1. 

Such a geometry may be generated by electromagentic fields only if these are considered in its average properties 
[p],pTp^ . Using the standard spatial average process we set m 



<E,> =0 
<B,> =0 



< E,Ej > 



< B,B, > 



< E,Bj > = 



-E' 



-B^9^J , 



Equations (h3^ - (17) then imply [m 



where 



< F^^F'', >= '^{E^ + B^)U^,U, + ^{E^ - 2B^)g^, , 



^^~ ds 



(13) 
(14) 

(15) 

(16) 
(17) 

(18) 

(19) 



Hence, for the classical Lagrangean (^ the average value of the energy-momentum tensor reduces to the form of a 
perfect fiuid 



where the density p and pressure p have the well known form 

P=\p- 



(20) 

(21) 
(22) 



In order to analyse the modifications implied by the use of the modified Lagrangean (j^), we assume that the 
dominant material content is a primordial plasma in which only the average value of the squared magnetic field B^ 
survives, i. e., we use Eqs. (|l^) - (|^ with E^ — ^. Then Equation (|^ still holds, but the energy density and 
pressure are now given by 

p=^B^l-2fiB'), (23) 

p^^B^l-10f,B')^^p~jf,B*. (24) 
/^From Eqs. (M), (ESf) and (p4) we see that the weak energy condition p > is obeyed if 

B < ^— = 6.0 X 10"''^ Tesla = 0.60 Gauss , (25) 

whereas the pressure will reach negative values only if 

B >2.7 X 10"^ Tesla = 0.27 Gauss . (26) 

Note, however, that expression (^) will hold only if the condition (|l^) is satisfied. Therefore, the conditions ( |lO|) and 
( p6| ) are not compatible. We shall return to this point later. 
The Einstein equations for the metric (|lj) read 

d A(t) AttG , 

- T^-^^(P + 3p), (27) 



a 



a^ K 8ttG A{t) 
from which the energy conservation equation can be written as 

P + 3-{p + p) = --^. (29) 

a 8ttG 

In the above equations we have allowed A to be time-dependent and the overdot means derivative with respect to the 
cosmic time t. 

Recent observations that have led to the so called age problem [13 - 15] have driven attention to cosmological 
models with a non-vanishing A. Estimates of the density parameter and some kinematical tests also point to the 
existence of an effective vacuum component [ [l6| , |l7| . Even more recently, some evidence for an accelerated cosmic 
expansion has arised from measurements involving type la supernova at high redshifts [HSlH. On the other hand, 
quantum field theory and inflationary cosmologies have pointed out to the possibility of treating A is a dynamical 
quantity [4,16^20 - 22] Several mechanisms have been identified as possible sources of fluctuating vacuum energy (see, 
for example, [Q and references therein). 

As it will be shown below, the field equations lead to a possible solution corresponding to a constant magnetic field 
plus a time-dependent cosmological term. In this case the time dependence of A is univocally determined. 

Replacing (pSh and (pi) in the Einstein equations (EtI) - (|2l) and in the energy conservation equation (p^), we get 

a A An G „o , „9s 

-^ B^il-GfiB'), (30) 

ad 3 

- + - = -^^2 (1-2^^2)+ (31) 

Qi Qj fj O 

B{l-4„B=)(B + 2£i,).-A,. ,3.) 

By solving any two of the above equations we find the complete cosmological solution for our model. 



III. GENERAL SOLUTION FOR A = AND A TIME-DEPENDENT MAGNETIC FIELD 

This is the case studied in R] where a particular solution for the scale factor was found. Wc will find the corre- 
sponding general solution and indicate how to recover this previous result. 

If B is time-dependent and A is a constant, Equation ( p^ ) can be easily integrated to give 

B{t)^Bo{^y, (33) 

where Bq is a constant of integration. In this paper, the subscript does not indicate the present day value of a 
quantity. Rather, it indicates the value of that quantity in an arbitrary time to, which will appear in the general 
solution for a{t) as a second constant of integration. It is this second constant to that was arbitrarily chosen in [nl. 
Thus Bo — B{to) and ao — a{to). (In Rcf. |l|] the authors normalized ao, so that oq = 1.) 

Note that (BSh holds for an arbitrary constant A. li A = K = 0, Eqs. ( pl| ) and (BSh allow us to find a{t). The 
general solution is 

a{t) = ao [4 al {t - tof -h 4 ao /3o (t - to) + l] ^'^ , (34) 

where we have defined 



UnG , , 

ao = V ^5— Bo , (35) 



(io = ^\-2^iBl. (36) 

In order to compare with the results of |1|] we recast (Q) in the form 

a{t)^ao{^aie + 4ao7ot + ^o)'/\ (37) 

with 

70 = /3o - 2aoio , (38) 

(5o EE \aoto [aoto - /3o) + 1 . (39) 

The linear term in t inside the parenthesis of ( p7| ) does not appear in the solution given by the authors of Ref . [Q . 
This is due to the fact that they have arbitrarily chosen the constant of integration to to be 



*°-2^-2B^V 47g ■ ^^°) 



^From the solution above we may determine the time behaviour of the magnetic field 

_Bo 
{4alt^ + 4ao7ot-f<5o)i/2 



m = .,„.2.. ,.„„., .^./. ■ (41) 



The energy density and pressure are obtained through the use of Eqs. (p3|), (p5), respectivelly. 
Note also that the Hubble parameter is 



a ao [2aoit-to) + Po 

n = 

Hence, 



a [4al{t - toy + iaoPoit - to) + 1] ' ^ ' 



Ho ^ mto) ^aoPo- Bo J tl^i^_ll^ . (43) 



(In |l|,^, the notation H is used to represent the magnetic field.) 



^From (^) we see that, for large t, we recover the usual solution for radiation dominated universes, a{t) oc i^/^. 
However, the more interesting feature of ( p7[ ) is that the quadratic function inside the parenthesis does not have real 
roots, being positive for any t. Therefore, the model is non-singular with a{t) reaching the minimum value 



ao (2 /i 52)1/4 (44) 



at 



+ 7o , /?o ...X 

tmin — —7^ — to — 7^ ■ (.40 j 

The universe thus obtained is a bouncing one: it begins arbitrarily large at i ^ t^^in, decreases until the minimum 
value (HI) at tmm and then begins to expand. 

A crucial feature of the above results is that the magnetic field at tmin is 

S(i„„0 = ^= = 6.0 X 10"^ Tcsla = 0.60 Gauss , (46) 

a value much greater than the minimal critical value given by Eq. (HQ) which limits the domain of validity of the 
first order correction given by (^). We are forced to conclude that the above solution does not apply all the way 
back to the point of minimum size of the universe. This point has not been acknowledged explicitly by the authors 
of reference ||l| , although their numerical calculations [g[^ based on the exact effective Lagrangean density (beyond 
the first order approximation) seems to confirm the non-singular behaviour of the model. We should remark that, in 
[Q, the requirement (^) is mentioned, but not the condition ([lO|), which is established only in equation (6.84) of [[lO[ . 
Note also that 

Pitrmn) = . (47) 



The times for which the constraint (10) is obeyed may be easily determined using ([411). We find that the solution 
(p7|) fulfill the above requirement for 



t < tjnin — t{l) and t > tjnin + t{l) , (48) 



where 



where TUpi = G^^/^ ~ 1.2211 x 10^^ GeV is the Planck mass. 

This seems to be a disappointing conclusion about the capability of the first order quantum correction (H) to 
describe the very early stages of the universe. In the next two sections we try to evade this problem by allowing for 
the existence of a non-vanishing cosmological term, as it has been increasingly indicated by recent theoretical and 
observational results U% . 

One could speculate that models with a time varying fine structure "constant" could perhaps evade the above 
mentioned problem. In fact, if a, and consequently /i, was grater at early times, the value of B{tmin) could reach 
values below the constraint (|l^) . However, in spite of the fact that models with oscillatory variation of the fundamental 
constants have been proposed |24U23], recent calculations using the Keck Telescope data seems to indicate a negative 
variation of a at high redshift (2 > 1): A a/a = -1.1 ± 0.4 x 10"^ ©l^l- 

Before going to this next step, let us mention that, if (5^ could describe the entire evolution of the universe in the 
distant past, it would imply the existence of an inflationary era (ii > 0) in the interval 

where 



t, ^ !!M . /^Z ^ 7.1 X 10^2 GeV-i = 0.046 s 

2 V 71" 



(51) 



Figure 1 shows the scale factor, the magnetic field, the energy density and the pressure as a function of time for a 
definite value of Bq. The time interval where the weak field approximation breaks down is also indicated. 



IV. GENERAL SOLUTION FOR A CONSTANT NON- VANISHING A AND A TIME-DEPENDENT 

MAGNETIC FIELD 

We will now extend the analysis made in ^ and investigate how the above conclusions are modified with the 
introduction of a constant cosmological term. 

For K = and A = constant ^ 0, substitution of ( |33| ) into (|3|) leads to the equation 



Z^ = 16 [XZ^ + al{Z~ 2^iB^)] , 



where we have defined 



"-U' ' 



(52) 



(53) 



A 



Equation (m) can be easily integrated to give 



a{t) = "0 ( ^ 



1/4 



C„e4^(*-*«) + ^e-4VA(t-to)_2, 



Co 



1/4 



where 



Co = al + 2X + 2^\{\ + al-2al^iBl) , 



(54) 



(55) 



(56) 



The Hubble parameter is 

H{t) = 2BoVA 



Do = al{al + 8\liBl) 



Co 



-1/2 



(57) 



(58) 



It is straightforward to see that the term inside the square brackets of ( p5|) is positive for all t and that the scale 
factor reaches its minimum value 



^min — ^0 



2X 



al + 8XfiB^ - ao 



1/4 



at 



1 fDa 

tmin — to ^ = In 



8VA V^o 



(59) 



(60) 



As in the previous case, the universe bounces at tmin and, if the solution would effcctivelly hold near tmin, an 
inflationary phase would take place for all values of t such that 



Co^x* - Sa'oCox^ + UDqx' - 8a^-f x + -^ > 



2^0 Dl 



where 



Nevertheless, the magnetic field at tmin is 

^\^7nin) 



„4VA(t-to) 



A 



27rG (,/l + %^-l 



1/2 



(61) 



(62) 



(63) 



^From the above expression we see that B{tmin) -^ —k= as A -^ and that Bitmin) ^ oo as A — > oo. Therefore 
this model suffers from the same problem as the previous one: the first order approxiamtion (p[) can not describe the 
dynamics all the way back to the point of minimum compression. 

^From the condition B{t) < Bcr, we find that the solution is valid for any time t such that 



t~U. 



In 



^2^ v'4 7rG(7rG + 2^A) 



where 



A = K + 2ttGBI.±^/K [A + 4 7rGB2,, (f-2/^B2J] 



(64) 



(65) 



The domain of validity of the solution will then depend on the value of A. 

Figure 2 and Figure 3 show the scale factor, the magnetic field, the energy density and the pressure as a function 
of time, for some values of A and Bq . 

V. GENERAL SOLUTION FOR A CONSTANT MAGNETIC FIELD AND A TIME-DEPENDENT A 



We will now turn to the case when the magnetic field does not vary with time. This solution of Eq. (p2) has not 
been analysed in Ref. R]. If A = 0, this leads to the well known models p = p = Oorp= —p — constant. However, if 
we treat A is a dynamical variable, then (p2) with 



B{t) = Bq = constant , 



leads to 



(66) 



A = 3^0- 



(67) 



and hence 



where 



and 



A{t) = Ao + 3Kq In 

' ao 



Ao = A(io) 



16 ttG o , o, 

Ko^ — i?2(l-4Aii?o') 



(68) 



(69) 



(70) 



Substituing this result into (|3^), we get for the scale factor 



a(t) — qq exp 



where 



{t - toY + Hoit- to) 



ao exp I -— t"^ + f3i t + /3o 



Ho^J^ + '-^Biil-2,Bi) 



Ao SttG 



■P , 



(71) 



(72) 



/3i 



Ko 



Ho 



(73) 



and 



Po = to[^t„-Ho] . (74) 



The Hubble parameter is 



and we have H = ior 



H{t) = ^it~to) + Ho 

= ^t + (3,, (75) 

tc = — ^=to--rr^. (76) 

At this point, the scale factor reaches the value 

a{t,)=aoe-"o/Ko ^ (77) 

whereas the cosmological term is 

A(tc) = Aq-3H^ = -At: G B^ {I - 2 fi B^) = -8tt G p . (78) 

The behaviour of the solution will depend on the sign of the constant Kq (for Kq = we get the de Sitter solution). 

For Kq > 0, that is, for Bq > 1/(2^/7*), the universe begins arbitrarily large as t — > — oo, contracts until reaches its 
minimum size at tc and begins to expands without limit from this point on. The solution is always acelerated (a > 0). 
Note that this case violates condition IfiVj). 

A much more interesting solution is the one corresponding to iiTg < (Bq < l/{2y/JI) w 4.2 x 10~^ Tesla = 
0.42 Gauss), since we can make Bq satisfy (|lO|) in this case. For this range of Bq, a{t) approaches zero asymptotically 
as < — » — CX3, reaches a maximum at tc and begins to contract, approaching zero as i — > +cxd. Moreover, a > for 



t<i.-y-|- and i>tc + \J-i-^- (79) 

It is worthy remarking that the time interval Ai^^v/); prior to t^ for which the solution is not inflationary depends 
on the value of Bq as 




This time interval reaches a minimum at Bq = l/(8/i) and goes to infinity as Bq —> or Bq -^ l/(4/i). Remembering 
that, for consistency, condition (nW must be obeyed, we see that this time interval may be as small as 3.2 s. Hence 
we have a cosmological model that could describe a universe that inflates through much of its history. 

Note that, as far as A{t) is concerned, tc is a point of minimum for both Kq > and Kq < 0. For the case of interest, 
{Kq < 0), the evolution proceeds in the following way (note that the weak energy condition (p > 0) is automatically 
guaranteed if Kq < 0). In the far past, the universe is arbitrally small, with an arbitrarilly large vacuum energy. It 
expands in such a way that during this expansion phase, the vacuum energy decays so as to keep p constant. The 
vacuum energy eventually becomes negative at the time 



2 Irr2 Ao 



reaching its minimum value 



^^^^^ + ]^V^°-f ' ^''^ 



Pvac - g^^ l»^j 



at tc, where the universe reaches its maximum size. Thereafter, the contraction phase takes place. The vacuum energy 
then grows, becoming positive at 



XoV ° 3 



Hl-^, (83) 



whereas p remains constant. The expected tendency of p growing adiabatically with the decreasing a is compensated 
by the increase in p^ac- The universe would end up arbitrarily small as i — > cxi with an arbitrarilly large vacuum 
energy. 

Figure 4 shows the scale factor and the cosmological term as a function of time for Xq > and some values of Bq 
and Ao, where we have further defined 

Ao ^ ^ . (84) 

Figure 5 shows the same quantities for A'o < 0. 

If one expects any such model to properly describe the evolution of the real universe, it would be advisable to take 
into account other matter fields, such as ultrarelativistic matter, scalar fields or dust. In |l||, it was demonstrated, 
for the case A = 0, B = B{t), that the presence of ultrarelativistic matter with an equation of state P{ur) = P{ur)/'^ 
would just amount for a reparametrization of the constants Bq and p. The effects of other matter fields in the models 
analysed in this section are currently under investigation. 

VI. SUMMARY AND CONCLUSIONS 

We have examined some consequences of considering first order quantum corrections to Maxwell elctrodynamics 
in zero curvature FRW universes. We have derived general analytical solutions under three different assumptions. 
When the cosmological term A is identically zero and the dynamics is driven by a time dependent magnetic field, we 
have obtained the general form of the solution found previously by Novello et al. [l| . We have found a new solution 
when A is a non- vanishing constant. In both cases, the universe is non-singular, bouncing at a critical time when it 
reaches its minimum size. However, we have shown that, near this critical time, the magnetic field increases beyond 
the value allowed by the use of the first order approximation to the effetcive Lagrangean density. The time range 
where this weak field approximation does apply was evaluated. On the other hand, numerical work, that takes into 



account higher order terms, indicates that the non-singular behaviour is preserved [p 23|. 

We have derived a third solution that describes a universe driven by a time dependent A that sustains a constant 
magnetic field. In this case, the time behaviour of the cosmological term is univocally determined and depends on 
the logarithym of the scale factor. For small enough values of the magnetic field strength (so that the first order 
corrections can be used safely), the universe begins arbitrarily small as i — > cxd, expands to a maximum size a„iax 
at tc and then contracts back to zero size as f ^ cx). The energy density associated with A is arbitrarily large at 
i — > — oo, reaches a minimum at tc and grows without limit as i — > c». It therefore decays during the expansion phase 
and increases during the contraction era, so as to keep the magnetic energy constant. The dynamics is inflationary 
during most of the cosmic history, except near tc- 

We are presently attempting to generalize our results for universes with non-zero curvature as well as analysing 
how our models would be modified by the presence of other matter fields. Another possibility is to study how the 
time evolution for the magnetic field, to be derived from (|32|), would be modified by imposing a definite decaying law 
for A. This could be chosen among the several forms presented in the literature (see, for example, H]). More accurate 
results could possibly be found by considering higher order terms in the effective Lagrangean. 
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FIG. 1. The upper panel shows the scale factor (solid line) and the magnetic field (dashed line) while the lower panel shows 
the energy density (solid line) and the pressure (dashed line) for the model with A = 0. The vertical bars indicate the time 
interval tmin — t{i) < t < tmin + i(i), during which the constraint (10) is not obeyed. Bo has been chosen such that ^/TJIBq — 0.2 
and Bo/B^ = 0.5. 
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FIG. 2. The upper panel shows the scale factor (solid line) and the magnetic field (dashed line) while the lower panel shows 
the energy density (solid line) and the pressure (dashed line) for the model with a constant non-vanishing A. The values for A 
and Bo are such that ^/X/ao = 2 x 10"'* and v^2/ii3o = 0.2. 
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FIG. 3. As in Figire 2 but for \/A/ao = 5 x lO"'^ and ^/2JIBo = 0.05. 




FIG. 4. The scale factor (solid line) and the cosmological term (dashed line) for the model with constant magnetic field, 
time-dependent A and Ko > {y^2flBo = 1). In the upper panel \/Ao/ao = 1 and the lower panel is for \/Ao/oo = 0.5. 
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FIG. 5. As in Figure 4 but for Ko < {-^flilBo — 0.1). In the upper panel \/\o/ao — 1 and the lower panel is for 
Ao/ao = 0.5. 
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